Classical limit of transport in quantum kicked maps 
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We investigate the behavior of weak localization, conductance fluctuations, and shot noise of a 
chaotic scatterer in the semiclassical limit. Time resolved numerical results, obtained by truncating 
the time-evolution of a kicked quantum map after a certain number of iterations, are compared to 
semiclassical theory. Considering how the appearance of quantum effects is delayed as a function of 
the Ehrenfest time gives a new method to compare theory and numerical simulations. We find that 
both weak localization and shot noise agree with semiclassical theory, which predicts exponential 
suppression with increasing Ehrenfest time. However, conductance fluctuations exhibit different 
behavior, with only a slight dependence on the Ehrenfest time. 

PACS numbers: 73.23.-b,05.45.Mt,05.45.Pq,73.20.Fz 
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I. INTRODUCTION 

According to Ehrenfest 's theorem, the expectation val- 
ues of the position and momentum of an electron obey 
classical equations of motion. As long as the wave- 
function of the electron is a wavepacket with minimal 
uncertainties in momentum and position, the expecta- 
tion values are a good description of the quantum state. 
However, wavepackets disperse, and the Ehrenfest theo- 
rem looses its relevance after a short time. In a cavity 
with point scatterers, which split electron wavepackets 
into partial waves after one scattering event, this time is 
simply the elastic mean free time. In a ballistic cavity 
with chaotic classical dynamics, this time is the so-called 
"Ehrenfest time" te, which depends on the Lyapunov ex- 
ponent A of the classical motion in the cavityi^ For times 
longer than te , a classical description no longer holds and 
the wave nature of the electrons becomes visible. 

The wave nature of electrons is the cause of some strik- 
ing effects that are absent in classical systems. For trans- 
port through cavities coupled to source and drain reser- 
voirs via point contacts, these effects are weak localiza- 
tion, universal conductance fluctuations, and shot noise. 3 
In the limit that transport through cavities is ergodic 
(dwell time in the cavity is much longer than the time 
of flight through the cavity), the signatures of quantum 
transport are 'universal', independent of the cavity size 
and shape, and of the fact whether electron motion in- 
side the cavity is ballistic and chaotic or diffusive, with 
repeated scattering off impurities with size smaller than 
the electron wavelength. Random matrix theory provides 
a unified theoretical description of weak localization, uni- 
versal conductance fluctuations, and shot noise in ballis- 
tic or diffusive cavities^ 

If the electron motion is diffusive, the dynamics is fully 
quantum mechanical already at times much shorter than 
the time r erg required for ergodic exploration of the cav- 
ity's phase space. For ballistic cavities this is true in 
most practical applications as well — the Ehrenfest time 
te usually does not exceed the time of flight through the 
cavity — but there is no fundamental reason why te al- 
ways has to be small. The case of large Ehrenfest times is 



of theoretical interest, as it is one of very few regimes in 
parameter space in which one can observe differences be- 
tween signatures of quantum transport in ballistic chaotic 
and diffusive cavities. 

The most prominent effects of a large Ehrenfest time 
are found if te is larger than the dwell time td in the 
cavity. If te 3> id, quantum transport is deterministic, 
and shot noise is suppressed m& The suppression of shot 
noise has been observed experimentally by varying the 
dwell time td of a chaotic cavity^ and numerically, using 
a chaotic map as a model for a chaotic cavity&2ii£ The 
effect of a large Ehrenfest time on weak localization was 
first addressed by Aleiner and Larking Their theory pre- 
dicts a suppression of weak localization oc exp(— te/td), 
if classical correlations are taken into account properly. 11 
The same suppression was found in an independent cal- 
culation by Adagideli. 12 Experimental observation of the 
suppression of weak localization at large Ehrenfest times 
has been reported for transport through antidot arrays^ 
No semiclassical theory for the Ehrenfest-time depen- 
dence of universal conductance fluctuations exists. How- 
ever, semiclassical theories for weak localization and uni- 
versal conductance fluctuations for the limit td 3> te are 
essentially equal^*i£*i£ii£ as are diagrammatic perturba- 
tion theories for the same phenomena in diffusive cavi- 
ties, supporting the expectation that the Ehrenfest-time 
dependencies of weak localization and universal conduc- 
tance fluctuations will be equal as wellii& 

Direct numerical simulation of the effect of a large 
Ehrenfest time on quantum transport through two- 
dimensional chaotic cavities has been problematic be- 
cause of the prohibitively high computational cost of the 
simulations. The reason is that te depends only loga- 
rithmically on the product of the electron wavenumber k 
and the cavity size L, 



te = A 1 In kL. 



(1) 



For two-dimensional cavities, system sizes at which te ^ 
td cannot be simulated with present-day algorithms and 
processor speeds. In order to circumvent this problem, 
Jacquod, Schomerus, and Beenakker proposed to replace 
the cavity by a quantum map^* The map is 'opened', 
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so that simulation of transport properties is possible. 
Although a map has a one-dimensional phase space, a 
chaotic map shares many characteristics of the chaotic 
motion in two-dimensional chaotic cavities^ The re- 
duced dimensionality of the map's phase space made 
numerical simulations with larger Ehrenfest times pos- 
sible. For an open version of the quantum kicked rota- 
tor map, numerical simulations were reported for shot 
noise*£u^£i weak localizationjiii2I and universal conduc- 
tance fluctuationsi2ii£i2£ Simulation results for shot noise 
were in good agreement with the predictions of the semi- 
classical theory, 6 » However, for conductance fluctuations, 
no dependence on te was found, despite the fact that 
Ehrenfest times larger than the dwell time were con- 
sidered^ Whereas early numerical simulations of weak 
localization showed no Ehrenfest-time dependence^ we 
showed that there is a systematic decrease of the weak lo- 
calization correction to the conductance upon increasing 
te, consistent with the semiclassical theory^i 

The main technical innovation that allowed us to de- 
tect a systematic decrease of the weak localization correc- 
tion upon increasing the Ehrenfest time is that we looked 
at time-resolved numerical simulations: The map's time 
evolution is truncated after a time to, and weak localiza- 
tion, conductance fluctuations, and shot noise are mon- 
itored as a function of t ~ This procedure has two ad- 
vantages. First, it allows the ensemble average over the 
quasienergy e to be done analytically. (See Sec. [n] for 
technical details.) This made it possible to consider sig- 
nificantly larger ensembles than considered previously. 
Second, monitoring quantum corrections as a function 
of the 'truncation time' to allows us to determine the 
minimal time after which quantum corrections can oc- 
cur. In the semiclassical theory, quantum interference 
requires a minimal wavepacket to be split and reunited, 
which takes a minimal time 2te- A schematic diagram 
drawing relevant semiclassical trajectories for weak lo- 
calization and conductance fluctuations is shown in Fig. 
n (The diag ram for conductance fluctuations is taken 
from Ref. [IT] and modified to contain the effect of a fi- 
nite Ehrenfest time.) Not being a quantum interference 
effect, shot noise only requires wavepackets to be split, 
which happens after a time te- Comparison of the time 
when quantum effects appear (the 'onset time') and the 
rate of suppression of quantum effects as the Ehrenfest 
time te is increased, thus provides a quantitative test of 
the semiclassical theory. With such a quantitative test of 
the semiclassical theory, accurate simulations performed 
at smaller Ehrenfest times can still be meaningful. 

In this paper, we present detailed results of such time- 
resolved numerical simulation for weak localization, con- 
ductance fluctuations, and shot noise. For all three quan- 
tum effects, we analyze their Ehrenfest-time dependence 
(without truncation of the time evolution) and the onset 
times (obtained from simulations with truncation of the 
time evolution). The simulation results are presented in 
Sec. m together with the predictions of random matrix 
theory for time-resolved transport through open quan- 
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FIG. 1: Schematic drawing of relevant trajectories for weak 
localization (a), conductance fluctuations (b), and shot noise 
(c). The letters "L" and "R" refer to the left and right con- 
tacts to the cavity, respectively. The trajectories shown in 
panel (a) are for the weak localization to the transmission 
T; the trajectories shown in panel (b) are for the covariance 
cov(7?l,^?r) of reflections from the left and right contacts. 
Both weak localization and conductance fluctuations require 
a minimal dwell time of 2te. Shot noise requires a minimal 
dwell time te only. 

turn maps. Whereas we confirm our earlier conclusion 
that the numerical simulations for weak localization show 
a suppression oc exp(— te/td), consistent with the semi- 
classical theory, we find that simulations for conductance 
fluctuations show a very small increase if te is increased, 
the effect being small enough to be consistent with the 
simulation data reported in the Iiteraturei2ii2i2£ The on- 
set times for conductance fluctuations are more than a 
factor two smaller than the onset times for weak local- 
ization, which is incompatible with the notion that con- 
ductance fluctuations arise as wavepackets are split and 
reunited. 

A second goal of this paper is to show how classical 
correlations are taken into account in the semiclassical 
theory of Ref. 0- The importance of classical correlations 
- no quantum diffraction takes place to or from clas- 
sical trajectories with classical dwell time shorter than 
the Ehrenfest time — was pointed out in the analysis 
of simulation data for the quantum kicked rotator, most 
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notably the simulations for shot noised The original ver- 
sion of the semiclassical theory^ which did not include 
these correlations, predicted a suppression of weak local- 
ization oc exp(— 2te/td). It is only after accounting for 
the classical correlations that the proper exponential de- 
cay oc exp(— te/td) is recovered. (Classical correlations 
are taken into account correctly in the semiclassical the- 
ory of shot-noise suppression at large Ehrenfest times A) 
In addition it is demonstrated that this semiclassical the- 
ory is unitary, i.e., that no probability is lost. Our dis- 
cussion of the semiclassical theory for weak localization 
can be found in Sec. IIIII We conclude in Sec. IIVI 

II. TIME RESOLVED TRANSPORT IN OPEN 
QUANTUM MAPS 

Despite their different dimensionality, chaotic quan- 
tum maps are believed to have the same phemonology as 
(closed) chaotic cavitiesj 24 ' 2 ^ 2 9i 2 7 provided one restricts 
attention to the "universal limit" of times longer than 
the "ergodic time" . (The "ergodic time" is the time re- 
quired to explore the phase space.) Open kicked quantum 
maps have been used as a model for transport through 
chaotic cavities* 1 ^ because maps allow to simulate sys- 
tems with many transport channels with relatively small 
computational effort. For short times, qualitative differ- 
ences between maps and cavities exist, however, because 
the time evolution in maps is discrete, whereas the time 
evolution in cavities is not. 

Quantum maps operate on a finite state vector ip of 
dimension M. In the context of the quantum kicked ro- 
tator these states are discrete quasimomentum states or, 
alternatively, positions on a lattice with periodic bound- 
ary conditions. In the language of a cavity, the elements 
of ip can be thought of as points on the cavity boundary. 
The time evolution of such maps is discrete, and given 
by the Floquet operator 

^{t + l)=^{t). (2) 

For the kicked rotator, this type of dynamics is relatively 
easy to simulate numerically, which is why it is used to 
numerically model quantum interference corrections at 
large Ehrenfest times. 

In order to study transport, the system has to be 
opened. Hereto two consecutive sets of and TVr ele- 
ments i/j are chosen that corresponds to two 'leads'^ The 
initial condition for a transport simulation corresponds to 
ip localized at one of the 'lead' points. Escape from the 
cavity is modeled by recording the amplitude of ip at the 
lead points at time t and setting ip to zero afterwards. 
Formally, this corresponds to the constructior&22*22iSI 

5(e) =P[1- e te TQ] ^ e ie TP T , (3) 

where S is the scattering matrix corresponding to the 
map, e is quasienergy, P is a (JVr + Nt,) x M matrix 
projecting on the sites corresponding to the left (L) and 



right (R) contacts, and Q = 1 — P T P. The conductance 
coefficients Gll, Grl, Glr, and Grr are defined as 

Grl = tr SClS^Cr., 

Glr, = tr SCrS'Cl 

G ll = tr SC h S^C h -N h , 

Grr = tr SC R S^C K -N R , (4) 

where Gr projects on the channels of the right lead, 
whereas Gl projects on the channels of the left lead. Uni- 
tarity implies Grl = G L r = -G L l = -Grr. The dwell 
time corresponding to the map is td = M/(Nl + Nr.). 

We note that, with this method of opening the map, 
the leads end abruptly. Such abrupt changes lead to 
diffraction effects, similar to diffraction from sharp cor- 
ners at the lead opening of a cavity. The semiclassical 
contributions of some diffracting orbits of this type were 
calculated for closed and open cavities, see, e.g., Refs. 
l3lU.32H.33l They were found to be an important contribu- 
tion to shot noise at small channel numbers N for rectan- 
gular cavities !?"^ However, to the best of our knowledge, 
the total contribution of such orbits to shot noise, weak 
localization, and conductance fluctuations in chaotic cav- 
ities at large N are unknown. Therefore, we compare 
our numerical results to the existing semiclassical theory, 
which neglects the effects of corner diffraction. 

In this section, we will study the transmission and re- 
flection coefficients as a function of a truncation time to- 
The truncation procedure involves writing Eq. J2J as a 
geometric series of a time-dependent scattering matrix 
S(t), 

CO 

S( £ ) = X> fet S(t), (5) 

i=l 

with 

S(t) = PITQ]^ 1 TP T . (6) 

Instead of taking the full geometric series ©, we truncate 
the series after time to, 

to 

S(t<t ;s) = £>-*£(<), ( 7 ) 

4=1 

and study transport properties as a function of to. Al- 
though such a truncation procedure does not represent 
a physical system, it allows us to get theoretical infor- 
mation on the times involved in quantum transport phe- 
nomena. In particular, we will verify at what truncation 
time weak localization and conductance fluctuations will 
first appear. 

The quantities of interest are the ensemble-averaged 
conductance coefficients Gll, Grl, Glr, and Grr, with 
S replaced by S(t < to), and the fluctuations of these 
coefficients. With the truncation procedure, S is no 
longer unitary, so that one does not necessarily have that 
Grl = G L r = -G L l = -Grr. Hence, we need to con- 
sider all four conductance coefficients separately. If T is 
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time-reversal symmetric, one still has Grl = Glr, and 
only three conductance coefficients need to be considered. 
For to 3> tq, unitarity is restored, and one conductance 
coefficient is sufficient. 

It is also of interest to study the to dependence of shot 
noise. This also can be done by replacing the scatter- 
ing matrix by S(t < to), as was done for weak local- 
ization and conductance fluctuations. The resulting to- 
dependent Fano factor is given bj£ 

= N 1 + m L tr [ SC ^C L (1 - SC K &C h )} . (8) 

The ensemble average is taken by first averaging over 
the quasienergy e, and then over various lead positions. 
The quasienergy average can be performed explicitly, see 
Eq. (7J above. For the fluctuations, we consider the vari- 
ances of the transport coefficients, as well as covariances 
of different coefficients. The variances are taken with re- 
spect to the quasienergy e, in order to ensure that the 
fluctuations are entirely of quantum originA 2 ^ The vari- 
ance or covariance for fluctuations with respect to e is 
then averaged over different lead positions. Again, hav- 
ing the explicit energy dependence at our disposal, 
the calculation of the variance or covariance with respect 
to variations of e can be performed explicitly. 



A. Random matrix theory of time-resolved 
transport 

As a reference for our numerical simulations, in which 
we take the Floquet operator J- of the quantum kicked ro- 
tator, we consider averages and fluctuations of the trans- 
port coefficients for the case that T is a random symmet- 
ric unitary matrix, taken from the circular orthogonal en- 
semble of random matrix theory. In the limit of large M , 
2Vl, and TVr, which is relevant for the semiclassical limit 
we consider throughout this paper, such averages can be 
calculated using the technique of Ref. EE We then find 

(G a0 (t < t )) = {G al3 (t < t o )) cl +SG a0 (t < t )+O(Ar 2 ), 

(9) 

where the indices a and (3 can be taken to be L and R. 
The first term (G a p(to)) c \ is the classical average, 



(G aP (t<t Q )) cl = a M P — -N a S a0 , (10) 

with x = 1 — 1/td, and 8G a0 (to) is the weak localization 
correction, 



SG aP (t < to) = M{l _ x) 5 aP - - — + 2 ) ■ 

For the covariances we find similarly in the limit of M, Nr, 3> 1, 

(cov e [G aP (t < to), G jS (t < to)}) = 2 ^ a ^ 7 ^ ^ + 2 3 (N y S 0S + N,S aS + N s S a , + N s S Pl ) TZ 2 

A " V; (SvySps + SasSpjKa + OiM- 1 ), (12) 



where 



(N R + N L + iy 



Tlx = — - 2x t °~ 1 ((to + l)x 2 + 6x - 1 + 1) - + x zto-2 

1 + X 1 + X 

2 1 \ 

■-(1 - x) 3 4 + (1 - x) 2 (l + x)t 2 + -i (l - x)(l\x 2 + 20x - 1) + 2x(x 2 + 5x + 2) J , 

x x^°~^ 

TZ 2 = + 2x t0 + (t + l^ 2 *"- 1 + (t - l)x 2t0 , 

1 + x 1 + x 

f 1 ™2to — 1\ 

= x ( —— x to ~ l + — . (13) 

\l + x 1 + x J K ' 



In Eq. (|12|l the symbol cov e denotes a covariance taken 
with respect to variations of the quasienergy only, 
whereas the brackets (...) denote the ensemble average 
over T . 



B. The open quantum kicked rotator 



In our numerical simulations, we take the Floquet op- 
erator T of the quantum kicked rotator map. The map is 
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described in detail by Tworzydlo et ai, see Ref. I2H The 
matrix elements of the Floquet operator are 

T mn = (XU^UUX) mn (14) 

with 

TJ — i\ f-l/2 2irimn/M 
v rnn — Ik/ e , 

Y _A „-i j ^ £ cos(27rm/A/+0) 



Here if is the so called stochasticity parameter that de- 
termines the classical dynamics of the map. The region 
K > 7.5 is associated with classically chaotic dynamics. 
The parameter <f> determines the precise quantization of 
the map and has no effect on the classical dynamics. In 
this model the size of the matrix M is even. In the sim- 
ulations we set Nr. = N-^ = N. 

We used the Floquet operator lfH|) to study conduc- 
tance fluctuations and shot noise. The Floquet operator 
(|14|l can also be used to study weak localization using 
the following argument: The average conductance {G a p) 
consists of a classical contribution and a quantum cor- 
rection, 

{G a p) = G a p. c \ + SG a f3- (15) 



In this model M is even, but not a multiple of 3. The 
kick potential is given by 

V{6) = K cos (7rg/2) cos 6 + \-K sin {irq/2) sin 29, (18) 

where q breaks the parity symmetry of the models The 
parameter 7 plays the role of a magnetic held and breaks 
time-reversal symmetry. By comparing the transport for 
7 7^ to that of 7 = it is possible to investigate the 
weak localization correction as a function of time. This 
is the method used in Ref. |2l[ 

The advantage of the second method is that it does not 
require the cancellation of classical conductances of two 
different systems. The disadvantage is that it involves 
the additional parameter 7, which itself also affects the 
map's classical dynamics. However, weak localization is 
affected on the scale 7 ~ 1/N, so that the effect on the 
classical dynamics can be expected to be small in the 
semiclassical limit N ^> 1. 

For both models, the Ehrenfest time is given by 

t e = \- 1 \iiN, (19) 
up to an A^-independent constant. The dwell time reads 

n, = i (20) 



The classical contribution G Q ^ iC i scales proportional to 
the channel number N, whereas the quantum correction 
has no iV-dependence (except for a possible weak depen- 
dence on N through the Ehrenfest time). The quantum 
correction can be extracted by comparing average con- 
ductances at N channels and 2N channels, 

SG aP « 25G a0 (N) - 5G a p(2N) 

= 2(G aP {N))-{G a(j {2N)). (16) 

In order to avoid a spurious contribution to 8G a p from 
classical conductance fluctuations, we made sure that the 
ensemble averages for (G a p(2N)) and (G a p(N)) were 
taken for precisely the same classical dynamics (same val- 
ues of K , same lead positions). This method was used 
previously in Ref. ITU 

Alternatively, one can study weak localization by con- 
sidering maps with and without time reversal symmetry. 
The kicked rotator map (|14|) has time-reversal symmetry. 
A simple extension of Eq. (|14fl that breaks time-reversal 
symmetry is the so-called three-kick modeljSi 

T mn = (xnY*UYnx) mn (17) 

where 

y _ X cos(27rm/M) 
1 ran — u rnn& ; 

X -A p -i^V(2nm/M) 

n„ m = M-VSe^/V^™-™) 2 . 



C. Numerical results 

The numerical algorithm used to iterate the maps 
(|l*4"|) and l(T7|l can be accelerated by using a fast Fourier 
transformi&SiiSSiS The results are computed as a func- 
tion of time and not of quasienergy. This allows to trun- 
cate the series (0 at the maximal time studied, and to 
calculate the e-averages explicitly. 

Weak localization. The value of the time-reversal sym- 
metry breaking parameter 7 at which weak localization 
corrections to the conductance are suppressed in the 
three-kick model (|17f) scales inversely proportional to the 
channel number A^^i Hence, in order to compare weak lo- 
calization corrections for different channel numbers, the 
average conductance is calculated as a function of the 
product 7 AT. Figure |21 shows the result of numerical sim- 
ulations of the ensemble-averaged conductance. The en- 
semble average is taken over 1000 samples, choosing K 
randomly in the interval 10 < K < 11.5 and varying the 
lead positions. The figure also shows a Lorentzian fit to 
the simulation data for small 7 AT. The dwell time is set 
to be td = 5 and the parity-symmetry breaking param- 
eter q = 0.2. Results for different A^ are offset vertically. 
When determining the magnitude of the weak localiza- 
tion correction, we restrict our attention to the range 
7^ < 1, for which the 7A^-dependence of (Glr) exhibit 
a pronounced dip. The width of this dip seems to be 
(almost) independent on system size, in agreement with 
the predictions from random matrix theory^i However, 
as can be seen from Fig. |5J the shape of this peak is not 
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well approximated by a Lorentzian at higher values of 7. 
For 7 AT > 1 the average conductance typically continues 
to increase with 77V, but at a much slower rate. 

Figure El shows the difference <5Glr(7) = (Glr(O) — 
Glr(7)} for 77V = 0.7, as a function of the cut-off time 
to. (These results are averaged over 20 000 different re- 
alizations.) The standard case (without truncation) cor- 
responds to the limit to — > 00. Upon increasing N, the 
figure shows a systematic decrease of <5Glr as well as a 
shift of the onset of weak localization to larger truncation 
times. 

To study these numerical result quantitatively, we ex- 
amine the dependence of the onset of weak localization 
and its magnitude on the number of channels. As an 
operational definition of the onset time t on , we define 
t on as that truncation time to for which #Glr(£o) — 
0.1<5Glr(^o ~ > 00). In principle, onset times could also 
have been defined using 5Gll- We prefer to use <5Glr, 
since the latter are less impacted by the discreteness of 
the map's time evolution. The left panel of Fig. 01 shows 
t on as a function of ./V for the simulation curves shown 
in Fig. |3| as well as for similar curves calculated for 
7./V = 0.3 (time-resolved data not shown) and for simu- 
lation data taken at dwell time td = 10. The right panel 
of Fig.^shows the dependence of |<5Glr| on N. Figure^ 
depicts the same data as figure^ but for stochasticity pa- 
rameter K taken uniformly in the interval 20 < K < 23. 

According to the semiclassical theory, the derivative 
dt on /dhiN = 2/A, where A is the Lyapunov expo- 
nent of the map. Similarly, semiclassics predicts that 
din |<5GLR|/dln TV = — 1/Atd- The classical map corre- 
sponding to Eq. (|17|l is described in Ref . 0, and its Lya- 
punov exponent is readily calculated using the method 
described in Ref. 39. Since the magnetic fields that 
break time reversal invariance are classically small, the 
Lyapunov exponents should be computed for 7 = 0. Nu- 
merically, we find A = 1.56 for K « 10 and A = 2.04 
for K ss 20. Lines with slopes corresponding these Lya- 
punov exponents are shown in Figs.0]and[Sl We conclude 
that the TV-dependence of |£Glr| is well described by the 
semiclassical theory. The onset times deviate somewhat 
from the expected slope, however. Since this deviation is 
more pronounced for the shorter dwell time we attribute 
it to fluctuations in the values of the Lyapunov exponent. 
Fluctuations of the Lyapunov exponent have a stronger 
effect on quantum corrections at short dwell times than 
at large dwell times. 2 

Although our data are obtained in the same way as in 
Ref.[2U our conclusions are markedly different. Whereas 
Ref. |2l] finds no evidence of a systematic Ehrenfest-time 
dependence of weak localization, we conclude that there 
is a systematic Ehrenfest-time dependence of weak local- 
ization and that the simulation results for weak local- 
ization are consistent with the semiclassical theory. We 
attribute the differences to a lack of accuracy in the sim- 
ulations of Ref.lU Indeed, on average the simulations of 
Ref. l2ll do show a slight decrease upon increasing N, but 
the statistical uncertainties are too large to rule whether 
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FIG. 2: Average conductance for the three-kick model l|17^ as 
a function of the time-reversal symmetry breaking parameter 
jN. Curves shown for N = 25, 50, 100 and 200 are offset 
vertically. The stochasticity parameter K « 10, q = 0.2, and 
the dwell time td = 5. Inset: Same for a larger range of 77V 
values. 
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FIG. 3: Time-resolved difference between ensemble-averaged 
conductance at 7./V = and 7./V = 0.7 for the three-kick 
model 1171 . The dwell time is td = 5, while K m 10 and 
q = 0.2 



the decrease is systematic or accidental. 

Figure shows time- resolved data for <5Glr for td = 
12.8, q = 0.3, jN — 0.2, and K taken uniformly in the 
interval 10 < K < 11.5. Whereas the simulation data 
are consistent with the semiclassical theory for N > 20 
(see the insets in Fig. |5J , the simulation results for 5Glr 
show significant deviations from the semiclassical theory 
for smaller channel numbers. This is no surprise, as the 
semiclassical theory is known to break down for small N. 
Despite the differences between the semiclassical theory 
and the magnitude of the simulated weak localization 
correction 5Glr for small N, the onset times appear to 
depend linearly on IniV down to the smallest channel 
numbers considered in the simulation (N = 5). 

The second method of computing weak localization 
uses the scaling of the conductance with the number of 
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FIG. 4: Left: Onset time t on obtained from the time-resolved 
difference Glr(7) — Glr(0). Right: Weak localization correc- 
tion SGlr, together with a fit oc exp(— te/td)- The stochas- 
ticity parameter K uniformly chosen between 10 and 11.5; the 
dwell times td is set at td = 5 and to = 10. For to = 5, data 
are shown for = 0.7 (circles), and 7 TV = 0.35 (squares). 
For to = 10, data are shown for = 0.3 (diamonds) and 
yN = 0.15 (triangles). 
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FIG. 5: Same as Figure^] but for K randomly chosen between 
20 and 23. 



FIG. 6: Time-resolved difference between ensemble-averaged 
conductance at yN = and yN = 0.2 for the three-kick 
model 1171 . for K w 10 and to = 12.8. Left inset: onset 
times. Right inset: weak localization correction. 
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FIG. 7: Weak localization correction to reflection and trans- 
mission coefficients, obtained using Eq. I|16fl . Data shown are 
for K « 10 and to = 5. The average was taken over an en- 
semble of 80 000 realizations, except for N = 400, as described 
in the text. 



channels, see Eq. (|16fl . Time- resolved results for reflec- 
tion and transmission coefficients are shown in Fig.[7]and 
|H1 with K taken randomly from the interval 10 < K < 
11.5 and dwell times = 5 and tq = 10, respectively. 
For td = 5, an average was taken over 80 000 samples, 
except for the data point at N — 400. Due to the nu- 
merical cost of the calculations for N = 400 only 40 000 
samples were considered for K ss 10 and 20 000 samples 
for K w 20. For td = 10, the average was taken over 
20 000 samples for K w 10 and over 40 000 samples for 
K w 20. In all cases the statistical error of the weak lo- 
calization correction was estimated to be ~ 10 -3 . Again, 
we compare the data in Fig. to lines with slopes which 
are fixed by the Lyapunov exponent of the corresponding 
classical map. For this model the Lyapunov exponents 
are 1.69 when K w 10 and 2.37 when K w 20. As for 
the three-kick model, the results are consistent with the 



semiclassical theory, see the right panel of Fig. |5| Fig- 
ure also contains results for K taken uniformly in the 
interval 20 < K < 23 (time- resolved data not shown). 

The time-resolved data shown in Figs. [7] and [S] show 
some structure at times smaller than the onset time. We 
cannot rule out that this structure, which appears to per- 
sist despite statistical averaging, affects our quantitative 
conclusions for the onset times and weak localization cor- 
rections. Since there is no such small-time structure in 
the weak localization data obtained by varying the field 
7, we conclude that the small-time feature in the data 
of Figs. and [S] must be a non-magnetic-field-dependent 
quantum correction (e.g., resulting from diffraction at 
the contacts). The feature disappears quickly upon in- 
creasing the stochasticity parameter K. 

The largest channel numbers in our simulations - 
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FIG. 8: Weak localization correction to reflection and trans- 
mission coefficients, obtained using Eq. 11611 . Data shown are 
for K « 10 and td = 10. The number of realizations used to 
calculate the ensemble average is 20 000. 
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FIG. 9: Left: Onset time t on obtained from the time-resolved 
weak localization correction SGlk- Right: Weak localization 
correction <5Glr, together with a fit oc exp(— te/td). Data are 
taken with stochasticity parameter K uniformly chosen be- 
tween 10 and 11.5 (circles) and between 20 and 23 (squares), 
and for dwell times td = 5 (solid) and td = 10 (dashed). 



N = 200 and N = 400 for the scaling method, with an 
average over 80 000 and 40 000 samples (for K ks 10), re- 
spectively — are smaller than the largest channel number 
used in the simulations Ref. |^ The observed decrease 
of the weak localization correction is statistically signif- 
icant and systematic, but small. We could not obtain 
sufficiently accurate simulation data for higher number 
of channels. Increasing N at fixed numerical cost implies 
that the number of realizations in the average has to scale 
oc N~ 2 , so that the statistical error scales cx N. Obtain- 
ing simulation data for N = 800 at the same numerical 
cost as our N = 200 data, would mean that only 5000 
samples can be averaged over, increasing the error by a 
factor 4. At that point the statistical error becomes com- 




FIG. 10: Time-resolved shot noise Fano factors for the open 
kicked rotator, with K = 10 and td = 5. Left inset: Onset 
times, determined from the location of the maximum of time- 
resolved Fano factor as a function of the truncation time to. 
Right inset: Fano factor compared to a fit tx exp(— te/td). 
Results for K « 10 (circles), K ~ 20 (squares), td = 5 (solid 
lines), and td = 10 (dashed lines) are depicted in the insets. 



parable to the expected incremental decrease of the weak 
localization data, and the simulation data loose statisti- 
cal significance. 

Shot noise. Figure lTUl shows time-resolved data for the 
shot noise Fano factor F, taken for the one-kick model 
with stochasticity parameter K sa 10 and dwell time 
td = 5. The statistical average was taken over 20 re- 
alizations only, which is sufficient as the Fano factor is 
self-averaging for large channel numbers. A more quanti- 
tative analysis of the data for K w 10 (circles) , K m 20 
(squares), td = 5 (solid lines), and td = 10 (dashed lines) 
is presented in the insets. Onset times, determined from 
the maximum of the F(to) graphs, are shown in the left 
inset of Fig. ^3 together with a linear fit oc te + const., 
using the Lyapunov exponents obtained from the cor- 
responding classical map. The time of the maximum 
was defined as onset time since it is problematic to ob- 
tain reliable onset times using the definition used pre- 
viously. However, numerically obtained onset times for 
shot noise are significantly impacted by the discreteness 
of the map's time evolution, and do not show a smooth 
linear increase with N over the range of channel num- 
bers investigated in our simulations. The right inset of 
the figure shows values of the Fano factors for the limit of 
large truncation times, together with fits oc exp(— te/td), 
with the Lyapunov exponents A obtained from the cor- 
responding classical map. We conclude that the rate of 
exponential suppression of shot noise upon increasing N 
is consistent with the rate of exponential suppression of 
weak localization. 

Conductance fluctuations. We performed time- 
resolved simulations for the variance of the reflection 
from the left contact Gll, the covariance of reflection 
Gll and transmission Glr, and the covariance of reflec- 
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FIG. 11: The truncation time-dependence of var Gll, 
cov (Gll, Glr), and cov (Gll, Grr) for the one-kick rotator. 
The dwell time is td = 5, while K is chosen uniformly be- 
tween 10 and 11.5. The corresponding time dependence given 
by random matrix theory is shown dashed. The short-time 
behavior is magnified in the inset. 



tions from the left and right contacts. We have computed 
the variance of the conductance coefficient Gll, as well as 
the covariances cov (Gll, Glr) and cov (Gll, Grr) as a 
function of the truncation time to- Unitarity of the scat- 
tering matrix implies that var Gll = —cov (Gll, Glr) = 
cov (Gll, Grr) if to - * 00 , but this equality does not 
hold for finite to- Our simulations were done for dwell 
times td = 5 and 10, and for stochasticity parameters 
K between 10 and 11.5 and between 20 and 23. Results 
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FIG. 12: Left: Onset time t OYL obtained from the time-resolved 
conductance covariance data. Right: Conductance covariance 
cov(Glr, Gll). Data are taken with stochasticity parameter 
K uniformly chosen between 10 and 11.5 (squares) and be- 
tween 20 and 23 (circles), and for dwell times to — 5 (solid) 
and td = 10 (dashed). 



for td — 5 and K 10, which are representative for all 
results obtained, are shown in Fig. 1111 The number of 
realizations used for the averages shown here is 80 000 
for td = 5, except for the TV = 800 data point, for which 
only 40 000 realizations were taken. The statistical error 
of the variances was estimated to be ~ 10 -4 . Onset times 
calculated from the time-resolved covariance of Gll and 
Glr are shown in Fig. 1121 together with conductance co- 
variances. The scatter of the data points in the right 
panel of Fig. ^] for the short dwell time td = 5 is larger 
than our statistical error and reproducable. It is proba- 
bly an artefact of the discrete time evolution of the map. 
(Note that similar scatter exists for the weak localization 
data, although, in that case, it is obscured by the large 
systematic decrease of the weak localization correction 
upon increasing TV.) 

As can be seen from the left panel of Fig.^J the onset 
time depend linearly on In TV: r on w 0.55 In TV + const, 
for K » 10 and r on w 0.38 In TV + const, for K w 20. 
However, the slopes are a factor 2.2 smaller than for the 
weak localization data, both for K « 10 and K w 20. 
Further, the conductance variance {i.e., the large-io limit 
of the data shown in Fig. lllfl shows a slight increase with 
increasing TV. The total increase is less than 10% over the 
range of channel numbers considered in our simulations 
and falls within the statist ical uncertainty of simulations 
reported in Refs. l9l ll8l22L where it was concluded that 
varG is independent of TV. Simulations of conductance 
fluctuations for the three-kick model slow a less than 10% 
decrease of var G for the same range of TV, as well as onset 
times that increase significantly slower with TV than the 
onset times for weak localization in the three-kick model 
(data not shown). 

Clearly, the simulation data for conductance fluctua- 
tions are qualitatively different from the simulation data 
for weak localization; they differ both with respect to the 
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TV-dependence of the onset times and the ^-dependence 
of the magnitude of the quantum corrections. This con- 
tradicts the notion that the same interference processes 
(diffusons and cooperons and their generalizations to bal- 
listic systems) underly both weak localization and con- 
ductance fluctuations, although one should note that 
there is no semiclassical theory of the Ehrenfest-time de- 
pendence of conductance fluctuations yet. 



III. QUANTUM TRANSPORT THROUGH 
CHAOTIC CAVITIES: SEMICLASSICAL 
THEORY 




C 



R 



A semiclassical theory for the Ehrenfest-time depen- 
dence of the weak localization correction in chaotic cav- 
ities was first formulated by Aleiner and Larkin. 2 That 
theory predicts that the weak localization correction is 
suppressed oc exp(— 2te/td). However, the analysis of 
Ref.0does not account for all classical correlations: it ne- 
glects the notion that no quantum diffraction takes place 
for electrons that spend less than a time te in the cavity. 
In this section we show that accounting for all classical 
correlations gives a weak localization correction that still 
depends exponentially on the Ehrenfest time, but with 
a different exponent: SG oc exp(— te/tjj). While a part 
of our calculations has already appeared elsewherejii this 
section includes a calculation of the full conductance ma- 
trix G a p . This allows to verify that unitarity is preserved. 
As the possible loss of unitarity is known to be a problem 
in semiclassical theories, it is important to demonstrate 
explicitly that the method used here preserves probabil- 
ity. 

The system under consideration is a ballistic cavity 
with chaotic classical dynamics. The cavity is two- 
dimensional, and it is coupled to two electron reservoirs 
via contacts of width c?l and c?r, see Fig. El Ergodicity 
of the electron motion inside the dot is ensured by the 
condition c?L)dR *C L. The cavity and the contacts are 
considered in the semiclassical limit di^k, d^fe 3> 1, where 
k is the electron wavenumber. Thus, one has the separa- 
tion of length scales 1/k <C d^, c?r <C L or, equivalently, 
the separation of time scales (vk) <C r org <C td, where 
v is the electron velocity, T crg ~ L/v the ergodic time, 
and Trj the dwell time. The latter can be expressed as 



-i _ {dh + d R )v 
T ° ~ ttA, 

where A is the area of the cavity. We also define 

dR „ d h 



Ft 



dn + dh 



(21) 



(22) 



as the probabilities that an electron at a random point 
in the cavity will leave through the right and left con- 
tacts, respectively. Electron motion inside the cavity is 
assumed to be fully phase coherent and electron-electron 
interactions are neglected. 



FIG. 13: Schematic picture of the ballistic chaotic cavity of 
size L. The cavity is attached to two leads, labeled "L" and 
"R", and with width c(l and dR, respectively. The contours 
Clj C-l, and Cr drawn in the leads are used for the calculation 
of transmission and reflection coefficients, see the text. 



The classical motion in the cavity is assumed to be fully 
chaotic (when closed), with a Lyapunov exponent A. This 
exponent measures the rate of divergence of close phase 
space points. Since this divergence is due to scattering 
from the curved walls of the cavity one can expect that 
A ~ l/r crg . We neglect fluctuations in the Lyapunov 
exponent and assume that its value is uniform. This is 
done to simplify the calculation as much as possible and 
should not affect the main physical conclusions as long 
as relevant times are much longer than T crg . 

The Aleiner-Larkin formalism is based on a semiclas- 
sical description of transport. Central objects in their 
theory are the non-oscillating parts of the product of an 
advanced and retarded Green functions. The "diffuson" 
is composed of combinations of orbits with themselves 
(known as the diagonal approximation in the semiclassi- 
cal community) while the "cooperon" is composed of or- 
bits and their time reversal counterparts. The equation 
of motion for the diffuson, neglecting quantum phenom- 
ena, is the Liouville equation, 



-iuj + Lx V°(uj; 1,2) =5(1,2) 



(23) 



where j = (</>.,-, Rj) = 1,2 denotes the phase space co- 
ordinates, limited to the energy shell: R is the electron 
position while the angle <p represents is the direction of 
the electron's velocity v. The operator L\ is the Li- 
ouville operator, taken with respect to the coordinates 
(01, Ri). For a hard wall cavity, L = v ■ d/dH, with ap- 
propriate boundary conditions at the walls. The symbol 
5(1, 2) denotes a delta function on the energy shell, that 
is, 5(1, 2) ee 27n5(Ri - R 2 )S(0i - <fa). 

Quantum effects are introduced by including a weak 
additional random quantum potential, treating it as a 
perturbation, and averaging over its realizations. 2 This 
random potential leads to two different effects, which 
are expected in a quantum theory. The first is a diffu- 
sion in phase space, which mimics the fact that quantum 
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dynamics cannot separate small phase space structures. 
The second effect is 'quantum switching' between close 
phase space trajectories, which results in weak localiza- 
tion. The strength 7 q of these effects is then set so that 
quantum corrections will appear after the correct time, 
the Ehrcnfest time te. Phase space diffusion modifies the 
evolution equation for the leading order diffuson asm 



-iuj + L\ - 7 q 



if 



P>; 1,2) = 5(1,2). 



(24) 



the strength 7 q of the phase space diffusion term is esti- 
mated to be£ 7 q ~ X 2 /kv. It is instructive to note that 



this term, although small, breaks the time reversal in- 
variance which is a property of both quantum mechanics 
and classical mechanics (in absence of a magnetic field). 
In the presence of quantum diffraction, the leading order 
cooperon, C°, also satisfies Eq. (|2U). 

In addition to the presence of the phase space diffusion 
term in the evolution equation 1)24(1 ■ quantum diffraction 
leads to an quantum correction to the diffuson: T> = 
V° + AT). This quantum correction is referred to as weak 
localization, because similar quantum corrections are the 
precursor of Anderson localization in disordered systems. 
The weak localization correction for the diffuson reads 2 



A2?(w;l,2) =V°(u;l,2 



C°(m;2,2) 
2ttKv 



C >^%°(.-,1,2) 



2-nfw 



d32? (w;l,3)£>V;3,2) 



2iu] — Lq 



•7q 



C°(^;3,3) 

2-kTiv 



(25) 



where j = (4>j + ir, Rj) denotes the time reversal of phase space point j and v = m/2iiTi 2 is the density of states per 
unit area. In terms of classical trajectories, the weak localization correction AT) corresponds to a configuration as 
shown in Fig. HJi: AT) is the correction to the product of a retarded and advanced Green function arising from the 
combination of a trajectory that intersects itself at a small angle and a trajectory that avoids intersecting itselfiS^ 
If the intersection angle is sufficiently small, such a pair of trajectories has strongly correlated actions, resulting in 
a net quantum correction to the diffuson even after ensemble averaging. The pairs of trajectories before and after 
the 'interference region' are represented by the leading-order diffusons T)° in Eq. I|25|l . whereas the trajectories in the 
closed loop — an orbit conjugate to its time-reversed counterpart — are represented by the leading order cooperon 
C°. 

In order to calculate the weak localization correction to the conductance, the dot's conductance coefficients Gll 
and Grl are expressed in terms of the diffuson T), 



Grl = 2ttTlvv 



dli / — — cos 



ci. 



-7T/2 



2tt 



n) I dh 



R — 7r/2 



-^cos(0r-0 2 )P(2,1), 
Ztt 



(26) 



and 



Gll = 2irhvv I dli I 
Jci 



2tt 



■ cos 



a I dl 2 j 



2tt 



cos(</. L - <^ 2 )2?(2, 1) - 1 



(27) 



Here Gl, Gr , and G L are smooth contours following a cross section of the contacts, see Fig. ^1 The contour G L 
is taken a small amount further away from the cavity than the contour Gl, see Fig. ^| The angles </>l and </>r 
denote the direction of the outward-pointing normal to Gl, Gr, whereas the angle L represents the direction of the 
inward-pointing normal to G L . Since we consider DC transport, the frequency to has been set to zero. We refer to 
the Appendix for a derivation. Similar expressions, but with an ambiguity in the precise definitions of the integrals 
over the lead cross sections, were derived by Takane and Nakamura^ 

The leading order "classical" transmission and reflection coefficients are found by substituting T>° for T>. Assuming 
that dynamics in the cavity is ergodic for all relevant time scales, one finds2*i& 



(28) 



where the escape probabilities to the left and right leads were defined in Eq. I|22|l above. 

The weak localization correction to the conductance coefficients Grl and Gll is obtained by substituting Eq. I|25|l , 
for the weak localization correction AT) of the diffuson, into Eqs. I|2t)|) and (|27|) . We first consider the weak localization 
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correction to Grl- The first two terms in (|25|l contain paths which leave the dot at the point of entry, and hence 
these do not contribute to the transmission. Hence 



6G KL = -v 2 



/ dh cos(^l-^i)/ dl 2 — 



0R. + 7T/2 



COS(0R 



x / d32> u (2,3)£> u (3, 1) 



C°(3,3). (29) 



r 



The phase space point 3 can be viewed as the center of 
the 'interference region' as depicted in Fig. QJl. The two 
diffusons and the cooperon head into opposite directions 
in phase space, so that they will sample different parts 
of phase space along the trajectory. Thus, one can treat 
them as statistically uncorrelated^ 

Both the cooperon and the product of diffusons in Eq. 
Pty exhibit non-trivial correlations since the returning 
path is strongly correlated with the outgoing one, cf. Fig. 
GJi. In order to take into account these correlations one 
should examine what is meant by the phase space points 
3 and 3 in Eq. (|2*9^l . By virtue of the phase-space diffu- 
sion term in the evolution equation for the leading order 
diffuson, effectively, the points 3 and 3 do not need to be 
exactly time reversed phase space points, since quantum 
mechanics allows for a finite phase space uncertainty. An 
intuitive way to calculate the correlations between 3 and 
3 was shown by Vavilov and Larkin*si Instead of using 
the diffusion in phase space, Vavilov and Larkin average 
over a range of initial phase space points near 3 and fi- 
nal phase space points near 3. They show that the effect 
of the phase space diffusion is equivalent to such averag- 
ing (up to logarithmic corrections) if the 'size' of phase 
space area near 3 and 3 scales as k~ x / 2 . [To be precise, 
the range of coordinates R in the average is ~ (L/k) 1 / 2 , 
whereas the angular range is ~ (Lk)- 1 / 2 ] We adopt the 
procedure of Vavilov and Larkin and consider an aver- 
age over phase space points 3' and 3" that are within a 
distance of order k^ 1 / 2 of the phase space point 3. We 
can then follow the classical orbits which start at the 
points 3', 3". The phase space distance between these 
trajectories will diverge exponentially due to the chaotic 
dynamics. Since the trajectories start from phase space 
points at a distance ~ k^ 1 / 2 , it will take a time te/2 to 
reach a phase space distance of order L. Once the phase 
space distance between the trajectories is large enough 
they can be considered as totally uncorrelated. At this 
point, ergodic dynamics can be assumed. The loss of cor- 
relations happens on a time scale of 1/A around the time 
Te/2. Since 1/A C te in the parameter regime where the 
Ehrenfest time affects quantum transport, we can view 
this process as effectively instantaneous. 

Let us first consider the cooperon in Eq. J2HJ). As a re- 
sult of the phase space correlations described above, one 
cannot close the orbit from 3' to 3" for short times. In 
order to make this more quantitative, we denote by tj the 
time it takes for the classical orbit starting at phase space 



point j to leave the cavity through one of the two open- 
ings. Then, if tg < t e /2 we find C°(3, 3) = C°(3",3') = 0. 
On the other hand, if > te/2 we can propagate the 
phase space point 3' for a time te/2 and reach phase 
space point 4. Similarly we can propagate 3" backward 
in time and find, toward the phase space point 5. This 
leads to 



C°(3",3"') 



if ig < r E /2, 

C°(5,4) if t 3 >t e /2. 



When the phase space points 3' and 3" are averaged 
over, the cooperon will have contributions from various 
phase space points 5,4, which are their distant past (or 
future). The points 4 and 5 can be taken to be uncorre- 
lated and sample the phase space with uniform probabil- 
ity. Since the phase space point 3 is eventually integrated 
over, there will be contributions from a sizable fraction of 
phase space. Thus, we approximate the average contribu- 
tion to the cooperon by replacing C°(5,4) by its average 
value td j A. This leads to 



(C°(3,3)> = ^3-t e /2), 



A 



(30) 



where 9{x) = 1 if x > and otherwise. 

We address the product of two diffusons in Eq. I|29|l by 
considering the integral 



J(3) = -v 2 



dh 



r<f>' h +ir/2 

dh / tt- C0S (^l 

R+V2 ^ 







2 1 ^cos(0-0 2 )P°(2,3)2?°(3,l). 

Cr J<Pr-7t/2 Z7r 



Again, the integration over phase space points 3 and 3 
is replaced by an average over phase space points 3' and 
3" within a distance of order (L/k) 1 / 2 from 3. Since the 
diffuson connects the phase space points 3' and 3" with 
points at two different contacts, the product of diffusons 
must be zero if t 3 < te/2 where, as before, t 3 is the time 
it takes for the orbit at phase space point 3 to leave the 
system. For larger £3 we may, again, assume ergodicity, 
so that we find, after averaging over 3' and 3", 

1{3) = P R P L 6{t 3 -T E /2). 

Combining results, we find 

TD 



SGrl = - p R p L-j / d36(t 3 - T E /2)L 3 9(t 3 - t e /2). 



(31) 
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The Liouville operator in Eq. (|31[1 measures the rate of 
flow of probability density out of the integration range 
of the phase space variable 3£ The boundary of the inte- 
gration range is composed of the lead phase space points 
propagated backward for a time £e/2. This leaves only 
a fraction exp(— ^r-) of the size of the boundary at the 
lead. To estimate the integral we assume that each out- 
going direction is equally likely to be in the system when 
propagated backward. However, only points which also 
have £g > te/2 will have a non- vanishing cooperon, lead- 
ing to an additional factor of exp(— gjjS^)- Collecting con- 
tributions from both leads, we find 

5G R L--P R P Le - TE / TD . (32) 



It is of interest to compute also the weak localization 
correction (JGll- This will allow to check that probabil- 
ity is conserved. The calculation is similar to that of the 
<5Grl, with two important differences. The first differ- 
ence is that the weak localization correction of reflection 
is composed from two parts. In addition to the third term 
in Eq. I[25[) . there is a contribution from the second term 
in Eq. . We write these two parts as SG^l and SG^l , 
respectively, and calculate them separately. Substitution 
of the third term in Eq. l(2"5|) into Eq. l(2"7|) gives 



= -v 2 



dh 



^+71-/2 



2tt 



-COS(<^l - 4>l) 



dh 



0L+7T/2 



l-tt/2 



2tt 



COS(</> L 



x / d3P°(2,3)P°(3,l) 



La 



1 d 2 



C°(3,3). (33) 



The product of the two diffusons can be assumed to be uncorrelated with the cooperon. The calculation of the cooperon 
proceeds as for the transmission calculation. However, the behavior of the product of the diffusons contributing to 
reflection differ from that of the diffusons contributing to transmission. To see that, again consider the integral 



1(3) 



dh 



#1 

2tt 



l) 



dh 



2vr 



cos(0 L -0 2 )2?°(2,3)P°(3,l) 



If h > te/2 both diffusons become uncorrelated before 
leaving the system and the average value of this integral 
is therefore P^. However, if £3 < te/2 both diffusons 
will exit the dot through the same lead, and one finds 
1 = Pl- Hence, 

fP L ift3<r E /2; 
\ P 2 if t~ 3 > te/2. 

This is the second difference between the calculations of 
<5Grl and <5Gll- We then find 

(2) 

The contribution SG^ is obtained by substituting the 
second term in the right hand side of Eq. I|25[) into Eq. 
()27J) . The lead integral over 2 can be calculated, resulting 
in 

6G$ = v f dh r L+ " 2 ^cos(^l-0i)C°(1, 1) 

= P L e- TE/2TD . (35) 

Combining Eqs. and one finds 5G LL = 50£l + 
(2) 

<5Gll = — <$Grl, as expected. 



The result (|32[) shows that, according to the semiclas- 
sical theory, the weak localization correction to transmis- 
sion is suppressed exponentially. However, the exponent 
we find is different from that of Ref. 0, where it is re- 
ported that <5Grl = — PL-FWexp(— 2te/td). The reason 
for the difference with Ref. is that our calculation obeys 
the classical correlations following from the separation of 
phase space into a 'classical' and 'quantum' part, corre- 
sponding to (classical) trajectories of length smaller or 
larger than te, respectively. Quantum diffraction does 
not involve the 'classical' part of phase space. Calculating 
the cooperon and the product of diffuson propagators as- 
suming ergodic dynamics in the 'quantum' part of phase 
space only increases the weak localization correction SG 
by a factor cxp(te/te) with respect to the calculation of 
Ref. |2j. The reason that SG remains exponentially small 

- but with exponent exp(— te/td), n °t cxp(— 2te/td) 

- is that, according to the semiclassical theory, weak 
localization requires a minimal path length of 2te«. The 
fraction of 'quantum' trajectories that remain inside the 
cavity during the time period 2te is exponentially small, 
oc exp(— te/td), hence the exponentially small weak lo- 
calization correction if te ^> td . 

Sofar we have considered the weak localization cor- 
rection to the transmission. What about other quan- 
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turn interference effects, such as the transmission fluctu- 
ations? Takane and Nakamura have extended the semi- 
classical theory of Aleiner and Larkin to the conduc- 
tance variance varG, but for the limit te — > only. 17 
They found varG = 2(PlPr) 2 , in agreement with pre- 
dictions from random matrix theory^ Just as the semi- 
classical theory for the weak localization correction <5Grl 
was simpler than the semiclassical theory for <5Gll (see 
above), the semiclassical theory of transmission fluctua- 
tions takes its simplest form if applied to the covariance 
of reflection from the right contact Grr and reflection 
from the left contact Gll- Of course, unitarity implies 
varG = cov (Gll, Grr). In Ref. 0, only one contribu- 
tion to the conductance covariance cov (Gll, Grr) was 
considered. This contribution corresponds to the four 
trajectories shown in Fig.^D. For these four trajectories, 
a minimal dwell time 2te is required: classical trajecto- 
ries originating from each of the openings need to diverge 
and reunite, each of which takes a time Tg. Hence, fol- 
lowing the same phase space arguments as for the weak 
localization correction, we anticipate that their contri- 
bution to the conductance fluctuations depends on the 
Ehrenfest times as 

[cov (Gll, Grr)] 1 / 2 = [varG] 1 / 2 oc e -W™ 

if T E >T D . (36) 

However, there may be contributions to the conductance 
fluctuations other than that of Fig. (This possibility 
can not be excluded on the basis of Ref.[l7].) For example, 
there may be trajectories that involve small-angle inter- 
sections of three trajectories at the same point 424^ While 
it is not expected that the inclusion of such trajectories 
undo the suppression of the conductance fluctuations in 
the limit te 3> td, they may have a non- negligible con- 
tribution for te ~ td- 

A final note about the results presented in this section: 
In the calculation performed in this section we have re- 
placed dynamical functions by their averages. For exam- 
ple, the probability that a wave packet will exit through 
a given lead is (almost) 1 or for short enough times, de- 
pending on the classical dynamics, not Pl or Pr. Thus, 
we have actually estimated the ensemble average of the 
weak localization correction. One can expect that if a 
system is close to the classical limit there will be more 
fluctuations in the values of diffusons and cooperons. 
This may lead to large deviations of the weak localiza- 
tion correction from that of systems with small Ehrenfest 
times, and to large classical conductance fluctuations. 22 
In order to describe these fluctuations quantitatively, one 
needs a sample-specific theory for the dynamics of the dif- 
fuson, which covers the times of order te exactly. While 
this is an important theoretical problem, it is beyond the 
scope of this paper, and should not affect the eventual ex- 
ponential suppression of quantum interference phenom- 
ena for large Ehrenfest times. 



IV. DISCUSSION 

In Sec. [HI we reported results of numerical simulations 
for weak localization, conductance fluctuations, and shot 
noise of the open quantum kicked rotator. The numerical 
simulations for weak localization and shot noise are con- 
sistent with an exponential suppression oc exp(— te/td), 
in quantitative agreement with the semiclassical theory 
of Sec. IIIII The numerical simulations for conductance 
fluctuations show a small increase with increasing Ehren- 
fest time, not inconsistent with previous simulation data 
reported in the literatureAiSiSS Our simulations also give 
information on the minimal time at which quantum ef- 
fects occur. For weak localization and shot noise, these 
times are 2te and te, respectively, consistent with semi- 
classical theory. For conductance fluctuations, the onset 
time is less than half the onset time of weak localization. 
(Note that it is possible that some contribution for weak 
localization in reflection may appear after a time of te-) 

In order to explain the numerical simulations of shot 
noise, weak localization, and universal conductance fluc- 
tuations, the authors of Refs. I3la ll8ll22l proposed a phe- 
nomenological alternative to the standard semiclassical 
theory, referred to as 'effective random matrix theory'. 
Guiding principle for the effective random matrix theory 
is that quantum diffraction takes place between trajec- 
tories with dwell times larger than te only; trajectories 
with dwell time shorter than te build scattering states 
with transmission (exponentially close to) or 1 and do 
not contribute to shot noise, weak localization, or univer- 
sal conductance fluctuations. The importance of 'quan- 
tum' trajectories is described by an effective number of 
'quantum channels' N q ~ 7Vexp(— te/td)- Following Sil- 
vestrov et al.^ it was then proposed that the scatter- 
ing of 'quantum channels' is described by random matrix 
theory as long as N q 3> 1. The condition N q 3> 1 is 
generically met, even if te ^> tq. 

The effective random matrix theory not only correctly 
predicts the suppression of the ensemble-averaged shot 
noise at large Ehrenfest times — shot noise is propor- 
tional to the number of channels N, which is replaced 
by N q = iVexp(— te/td) in the effective random ma- 
trix theory — , it also describes sample-specific deviations 
from the ensemble average that arise from the classical 
dynamics. 8 The effective random matrix also has had re- 
markable success explaining other observables that are 
proportional to the channel number N, such as the den- 
sity of transmission eigenvalues^ as well as the density of 
states in a chaotic cavity coupled to a superconductor. 44 
On the other hand, quantum-interference effects, such as 
weak localization and conductance fluctuations are in- 
dependent of N, and, hence, are predicted to be inde- 
pendent of the Ehrenfest time. Thus, for weak local- 
ization, the effective random matrix theory differs from 
the semiclassical theory, which predicts a suppression 
oc exp(-T E /r D ). 

The effective random matrix theory and the semiclas- 
sical description of transport not only disagree regarding 
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the magnitude of the weak localization correction to the 
conductance, they also disagree with regard to the min- 
imal time required for quantum interference effects to 
occur. In the semiclassical theory, quantum interference 
requires a minimal wavepacket to be split and reunited, 
which takes a minimal time 2te- This is in contrast to the 
effective random matrix theory, where quantum interfer- 
ence is fully established already after a time Te. Inter- 
estingly, there is no difference between semiclassics and 
effective random matrix theory for shot noise: Not being 
a quantum interference effect, shot noise only requires 
wavepackets to be split, which happens after a time te 
in both theories. 

These two differences between the semiclassical theory 
for weak localization and the effective random matrix 
theory are not unrelated. Both theories are consistent 
with a fully classical description of electron dynamics for 
the first time interval of length te after the electron has 
entered the cavity. Differences between the two theories 
appear for electrons that escape from the cavity in the 
second time interval of length te (i.e., for times between 
te and 2te). In the semiclassical theory, electrons that 
escape during this interval behave quantum mechanically 
but do not contribute to weak localization. As explained 
in Sec. 11111 the escape of electrons between te and 2te is 
responsible for the suppression of weak localization with 
exponent exp(— te/td) in the semiclassical theory. In the 
effective random matrix theory weak localization sets in 
as soon as the classical description fails, one Ehrenfest 
time after the electron enters the cavity; there is no time 
interval in which electron dynamics is neither fully clas- 
sical nor described by random matrix theory. 

A 'microscopic theory' supporting the effective random 
matrix theory and its prediction of Ehrenfcst-time inde- 
pendent weak localization and conductance fluctuations 
appeared recently^ This theory is based on the assump- 
tion that electron dynamics in the quantum trajectories 
(i.e., trajectories with dwell times larger than te) is fully 
ergodic, with fully established quantum interference cor- 
rections. Such an assumption violates the semiclassical 
picture of weak localization, in which the quantum cor- 
rection only appears after a time 2te- Hence, a theory 
which tries to justify the effective random matrix the- 
ory must provide an entirely new semiclassical model for 
weak localization, in which the quantum correction ap- 
pears after one Ehrenfest time only. 

Although, at present, there is no theory that explains 
all simulation data, we can conclude that our time- 
resolved simulation data for the quantum interference 
corrections to the conductance are not consistent with 
the effective random matrix theory. First, because the 
effective random matrix theory predicts onset time te for 
weak localization, as well as an Ehrenfcst-time indepen- 
dent weak localization correction, both of which are ruled 
out by our simulations. Second, because simulation re- 
sults for weak localization and conductance fluctuations 
are qualitatively different, whereas the 'effective random 
matrix theory' predicts equal onset times and Ehrenfest- 



time dependences for weak localization and conductance 
fluctuations. 

Of course, the question why numerical simulations for 
weak localization and conductance fluctuations in the 
open quantum kicked rotator are qualitatively different 
remains. Clearly, this question can not receive a final 
answer as long as there is no semiclassical theory for the 
Ehrenfest-time dependence of conductance fluctuations. 
Given the difficulty of obtaining accurate numerical re- 
sults in the regime te 3> td , such a semiclassical theory 
needs to include the parameter range te ~ td if a valid 
comparison with the numerical simulations is to be made. 
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APPENDIX A: TRANSMISSION AND 
REFLECTION IN TERMS OF DIFFUSONS 

In this Appendix the transmission and reflection 
through a quantum dot are derived in terms of lead in- 
tegrals of a diffuson. Similar expressions were derived 
previously by Takane and Nakamura. 16 However, their 
expressions are ambiguous with respect to the exact lo- 
cations of the cross sections in the leads. In our approach, 
all ambiguities are resolved. 

To calculate the total transmission and reflection in 
terms of diffusons, it is useful to obtain expressions 
for scattering matrix elements using the retarded Green 
function of the cavity. Exact expressions of this type were 
derived in Refs. The leads have a uniform cross 

section, so that the lead wavefunction can be decomposed 
into free waves along the lead (denoted by x coordinates) 
and a basis of transverse wavefunctions Xa(y)- The same 
decomposition can also be used for the retarded Green 
function (at the Fermi energy), and one defines 

g+(ri,r 2 ) =^2g+ n {xi,x 2 )Xm{yi)Xn(y2)> ( A1 ) 

mn 

for coordinates ri and r 2 in the leads. We use the con- 
vention that different leads are assigned different modes, 
so that Eq. 1|A1[) remains meaningful if ri and r 2 are in 
different leads. 

Using the asymptotics of Q + and of the scattering 
states, the Green function in the leads can be written 
in terms of scattering matrix elements S mn . One finds 46 



G mn (xi,x 2 ) = 




(A2) 
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for x\ > xi with x\ % in the left lead, and to 



Q^ n {xi 1 x 2 ) 



n „ik m xi 



-ik n X2 



(A3) 



for x\ in the right lead while x 2 is in the left lead. (In 
the following we use units where % = 1.) For the follow- 
ing discussion it is important to emphasize the condition 
x\ > .x 2 , which appears when both x\ and X2 are in 
the same lead. One of the steps in the derivation of Eq. 
(|A2(I involves an integral over a finite domain (bounded 
by the cross-section x — a; 2 ), where the integrand is pro- 
portional to S(ri — r 2 ). This integral will give different 
results depending whether ri is in the domain of inte- 
gration (xi > x 2 ) or not (x\ < ir 2 ). This leads to the 
inequality in Eq. 1A2(1 . Later, when we write in 
terms of cross section integrals, two different cross sec- 
tions, at x\ and X2, are obtained. 

To compute the conductance coefficients Gll and Glr, 
the absolute value squared of the matrix elements is 
needed. It can be expressed in terms of Green function 
with the help of 



(— 


-—) 


(— 


-—) 


\dxi 


8x4 J 


\dx 2 


dx 3 J 



XI— X4 

x 2 =x 3 



4m 2 1 5, 
4m 2 (6, 



mn\ 2 if Xi £ R, X 2 £ L 

mn ~ \&mn\j if x l t x 2 S L, X\ > X 2 ■ 



(A4) 



This expressions differ from those in Ref. in that the 
integrals for the reflection are along separate cross sec- 
tions, with the cross section Cl closer to the cavity than 

On each of the cross sections there are two points which 
are almost identified: ri ~ and r2 ~ r3. Therefore, it 
is natural the use the following Fourier transforr^ ' ' ' 



G + (r 1 ,r 2 )g (r 3 ,r 4 



dpi 



dp 



2 gipi-Cri-r*) 



(2tt) 2 J (2tt) 2 

x e ^-(r3-r 2 )^ (pijRl . p2;R2)j (A8) 



where Ri = (r x +r 4 )/2 and R 2 = (r 2 +r 3 )/2. Note that 
Ri and R 2 are on the cross sections of the lead. Energy 
conservation can be used to simplify this probability den- 
sity further. This is done by defining 



The conductance coefficients Gll and Glr are then cal- 
culated using Eq. J3J. One can use Eq. (|A4|) to simplify 
the summations, since the lead modes enter into G^n y i a 
cross section integrals over the leads. One then uses 



^ 23 (p 1 ,R 1 ;p 2 ,R 2 ) = — S[E. 



2m J V 2m 



xZ>( ni ,R i; n 2 ,R 2 ). (A9) 



E 

n=l 



x* n {yi)xn{yi) ^ S\ F (vi - yi)- 



(A5) 



Replacing this approximate finite width delta function by 
an exact one leads, after a straightforward calculation, to 



Grl 



1 



4m 2 
d 
dx 2 



dyi / dy 2 



d d 
dxi 8x4 



d 

dx 3 



g + (ri,r 2 )[g + (ri,r 3 )Y 



Gll — — 



1 



4to 2 



dyi 



thi> , _d d_ 

dxi 8x4 



ri=r 4 
f2=r3 



(A6) 



ri=r 4 

r2=T3 



(A7) 



Substitution of l|A8jl and (|A9|) in (|A6|) will lead to the 
expressions for Grl and Gll in terms of the diffuson. 
The calculation is straightforward. Some care is needed 
only when one determines the boundary conditions. For 
Grl , only diffusions that go into the cavity can get from 
the left lead into the right lead. These diffusons can arrive 
to the right lead only if they originate from the cavity. 
These considerations result in the Theta functions in Eq. 
(|26|l . The calculation for Gll is very similar. The only 
difference is that both incoming and outgoing diffusons 
will cross the inner (Cl) cross section. However, all the 
incoming probability flux from C L must cross Cl- This 
incoming flux just cancels the factor in i|A7|l . The 
remaining contribution results in Eq. (|27|l . 
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